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HORIZONTAL  RAY  CURVATURE  EFFECTS  IN  BASINS,  TROUGHS, 
AND  NEAR  SEAMOUNTS  BY  USE  OF  RAY  INVARIANTS 


INTRODUCTION 

Long-range  underwater  sound  propagation  is  often  influenced  by  regions  of  repeated 
surface  and  bottom  reflection.  Under  these  conditions,  a component  of  bottom  slope 
consistently  to  one  side  of  the  ray  path  may  cause  considerable  horizontal  curvature.  For 
instance,  it  is  easily  shown  (by  considering  source  images)  that  the  horizontal  projection 
of  the  ray  path  above  a plane  sloping  bottom  is  a hyperbola.  Similarly,  in  a basin,  one 
mi^t  expect  to  find  ray  paths  that  loop  around  the  center  to  form  a rosette-like  shape. 

It  is  important  to  note  that  the  curvature  effect  applies  to  rays  of  small  as  well  as  large 
grazing  angles.  Despite  the  laj^er  number  of  reflections  requh^d  at  low  grazing  angles  to 
produce  a given  change  in  heading,  the  consequent  reduction  in  reflection  coefficient  may 
allow  propagation  losses  that  are  not  prohibitively  large.  This  phenomenon  was  first  dealt 
with  by  Weston  [1,2] , and  some  of  the  effects  in  the  vertical  plane  were  investigated 
later  by  Milder  [3] . 

Field  recordings  of  remote,  powerful  pulsed  sources  (sometimes  the  equivalent  of  a 
ton  or  more  of  explosives)  [4-10]  often  contain  discrete  and  continuous  arrivals  stretch- 
ing over  periods  on  the  order  of  10  min,  depending  on  the  size  of  the  basin.  This  type 
of  echo  is  often  interpreted  in  a phenomenological  way  as  a single  reflection  or  scatter 
from  one  or  more  seamounts  or  basin  walls  [4-6] . Whether  the  ray  deflection  is  really 
caused  by  multiple  surface  and  bottom  reflections  or  by  a single  scatter  remains  to  be 
proved,  but  whatever  the  cause,  as  long  as  the  seabed  is  not  horizontal,  after  bottom 
interaction  the  ray  path  cannot  be  considered  to  remain  in  the  same  vertical  plane. 

This  report  attempts  to  find  a detailed  theoretical  explanation  for  the  long  ray  paths 
encountered  in  basins  by  using  ray  invariants  to  calculate  the  horizontal  projection  of  the 
multiply-reflected  ray.  The  treatment  differs  from  that  of  Weston  [1,2]  in  choosing,  as  a 
starting  point,  a simpler  sound  velocity  structure,  but  more  complex  bottom  topognq)hy. 


I.  THREE-DIMENSIONAL  RAY  PATHS 

To  construct  a three-dimensional  ray  path,  it  is  first  necessary  to  calculate  the  change 
in  ray  an^e  after  a single  reflection. 

Figure  1 shows  a ray  striking  a bottom  that  is  assumed  to  slope  down  at  angle  7 in 
the  y direction;  0 , 0 are  incident  ray  elevation  and  azimuth  angles  with  a prime  added 
for  reflected  rays;  and  <!>  is  the  local  angle  of  incidence.  We  wish  to  find  expressions  for 
6'  and  4>'  in  terms  of  the  other  variables,  and  these  may  be  shown,  by  straightforward 
spherical  trigonometry,  to  be 
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Fig.  1 — The  geoqietry  for  a sin^e 
reflection 


sin  fl ' = cos  27  sin  0 + sin  27  cos '6  cos  0 


and 


tan  0'  = 


sin  0 cos  B 

(sin  0 sin  27  - cos  6 cos  27  cos  0) 


(1) 

(2) 


The  change  in  heading  of  the  ray  after  reflection  o,  (where  a = tt  - 0 - 0')  is  given 
from  Eq.  (2)  by 


^ tan  0 (sin  sin  27  - cos  27  cos  6 cos  0)  + sin  0 cos  d 
tan  0 (sin  0 cos  0 + cos  27  cos  6 cos  0)  - sin  0 sin  27 

and  this,  as  expected,  is  zero  if  the  bottom  slope  7 is  zero.  Another  useful  equation  that 
is  independent  of  7 may  be  derived  from  Fig.  1,  and  this  is 


cos  0'  sin  0'  = cos  0 sin  0. 


(4) 


Consider,  now,  a ray  in  a uniform  medium  bounded  by  a plane  upper  surface  and  a 
bottom  of  arbitrary  form  but  low  slopes.  At  the  ith  bottom  reflecting  point  we  denote 
local  slope,  incident  elevation  and  azimuth  ancles  by  7,-,  0j,  and  0,-,  respectively.  Since 
the  upper  surface  is  horizontal,  the  elevation  angle  of  the  ray  reflected  from  the  tth  point 
d'j  is  the  same  as  the  elevation  of  the  incident  ray  at  the  (i  + l)th  point  0^^;  however,  a 
similar  relation  is  not  generally  true  tor  the  azimuth  angles  0,'  and  because  the  direc- 
tion of  bottom  slope  may  change  from  point  to  point  (Fig.  2).  After  the  first  reflection, 
we  have  from  Eq.  (1) 


sin  01  = sin  0q  + 279  cos  Oq  cos  <I>q,  (5) 

which  for  later  convenience  we  write  as 

sin  01  = sin  0q(1  + 27o  cot  00  cos  0o).  (6) 
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Fig.  2 — Horizontal  projection  of  a ray  path 
ihowing  bottom  reflection  points  fi,  surface  re- 
flection points  S,  and  the  changes  in  bottom 
slope  directions  y in  between 


After  the  nth  reflection,  repeated  use  of  Eq.  (6)  gives 


n-\ 

sin  d„  = sin  0q  (1  + 27,-  cot  0,-  cos  0,). 


t 


Taking  logarithms  of  both  sides  yields 

n-l 

Iog(sin0„)  = log(sin0o)+  log  (1  27,  cot  0,  cos  0,).  (8) 

i-0 

Provided  27,-  cot  0,-  cos  0,-  is  much  less  than  unity,  the  term  log  (1  + 27,-  cot  0,-  cos  0,) 
may  be  replaced  by  27,-  cot  0,-  cos  0,-.  For  this  to  be  true,  the  rays  must  not  be  too  close 
to  the  horizontal,  although  this  restriction  is  not  severe  since  7,-  is  already  very  small. 
Rou^y  speaking,  ray  elevation  angles  must  be  greater  than  twice  the  local  bottom  slope 
so  that  rays  cannot  be  reflected  downward  from  the  bottom.  Equation  (8)  is  now 

n-l 

log  (sin  0„)  = log  (sin  0q) -(■  27,’ cot  0,- cos  0,-,  (9) 

i-0 

and  the  summation  may  be  expanded  in  terms  of  the  water  depth.  The  horizontal  spac- 
ing between  bottom  reflections  (Pig.  2)  is  As,  = 2Hj  cot  0,-,  and  the  component  of  As,  in 
local  downslope  direction  y,'  is 

Ay,-  = As,-  cos  0,-  = 2W,-  cot  0,-  cos  0,.  (10) 

Now,  since  7,-  = - dHjIdyj,  the  summation  in  Eq.  (9)  may  be  replaced  by 
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n-1 


E 

i-O 


Ay, 


and,  provided  the  bottom  slope  changes  smoothly  between  one  reflection  point  and  the 
next  so  that  the  term  in  brackets  varies  smoothly,  we  may  replace  the  summation  sign  by 
an  integral  sign.  Then  we  have 


log(sin0„)  = 


log  (sin0o) 


4,  H 


which  leads  to 


sin  dn  = sin  do 


Hr, 


(11) 


(12) 


Thus  H sin  d is  invariant  for  a particular  ray,  and  this  is  a special  case  of  Weston’s 
invariant  [1] 


For  a given  initial  elevation  angle  dg,  the  elevation  angle  at  another  site  depends  only 
on  the  ratio  of  the  water  depths  at  that  site  and  at  the  source  (regardless  of  the  coordinate 
system  or  functional  form  at  the  seabed).  On  a map  of  bottom  topogn^hy,  contours  of 
water  depth  are  also  contours  of  ray  elevation  angle. 

This  result  could  have  been  derived  much  more  directly  from  Eq.  (1)  except  that  the 
approximation  used  between  Eqs.  (8)  and  (9)  would  not  have  appeared  so  clearly. 

Another  result  may  be  derived  from  Eq.  (4);  however,  this  time  the  validity  does  not 
depend  on  small  bottom  slopes  although  a usable  form  depends  on  the  shape  of  the 
bottom. 

Using  Eq.  (4)  repeatedly  for  a trough  of  arbitrary  but  constant  cross  section,  so  that 
the  y,  axes  all  point  in  the  same  direction  and  and  0,'  = jt  - we  have 


cos  sin  0,+j  = cos  0,  sin  0,; 


(13) 


and,  after  n reflections. 


cos  0„  sin  ^ = cos  0o  sin  (pQ. 


Thus  for  a trough  of  fixed  cross  section,  cos  6 sin  (p  is  invariant. 


(14) 


NRL  REPORT  8144 

For  the  more  general  case,  as  in  Fig.  2,  where  maximum  bottom  slope  direction 
is  a function  of  position,  we  have  and  ip’j  = ir  - - A4>,+i , where  A4>,+i  is 

the  change  in  maximum  bottom  slope  direction  encountered  between  the  ith  and  (i  + l)th 
bottom  reflection  (assumed  to  be  a small  angle). 

At  the  ith  reflection. 


cos  6'.  sin  0!  = cos  0,-  sin  0,-;  (15) 

and  at  the  (i  + l)th  reflection, 

cos  sin  0/^j  = cos  d._^^  sin  (tt  - 0,>i  - A4>,+i ) = cos  0,+i  sin  0,+i  (1  + cot  0,>i  A<l>,+i ). 

(16) 

Using  Eq.  (16)  repeatedly  and  taking  logarithms  of  both  sides  results  in 

n 

log  (cos  dl,  sin  0,^)  = log  (cos  0o  sin  0^)  + ^ log  (1  + cot  0,  A*!*,).  (17) 

1=0 


In  the  last  term,  cot  0^  may  be  written  as  (Ay^/AjCj),  where  Ay;  and  Axj  are  the  y and  x 
components  of  the  distances  between  the  ith  and  (i  - l)th  reflection  point.  Provided  the 
ray  is  not  close  to  the  intersection  of  the  y,  and  y,+i  axes,  which  is  the  local  center  of 
curvature  of  the  depth  contours,  A<I>,/Ax,-  is  a good  approximation  for  the  contour 
curvature,  and  the  last  term  in  the  summation  is  much  less  than  unity.  If  we  use  the 
small  argument  approximation  for  log  (1  + X),  we  have 

log  (cos  sin  <!>'„)  = log  (cos  Oq  sin  0o)  + 2^  - Ay,,  (18) 

where  p is  the  local  radius  of  curvature  of  the  depth  contours. 

If  the  curvature  of  the  depth  contours  changes  smoothly  from  reflection  point  to 
reflection  point,  the  summation  sign  can  be  replaced  by  an  integral  sign: 

log  (cos  0,^  sin  0j,)  = log  (cos  00  sin  00 ) + . (19) 

If  the  change  in  position  of  the  center  of  curvature  is  always  much  smaller  than  Ay,, 
then  dy  may  be  replaced  by  -dp  and  integration  performed  to  give 

cos  91,  sin  0j,  p„  = cos  6q  sin  0o  Pq.  (20) 

However,  this  last  condition  is  somewhat  restrictive,  and  the  only  type  of  basin  for  which 
Eq.  (20)  is  unconditionally  true  is  that  with  either  parallel,  straight  depth  contours  or 
concentric,  circular  contours.  The  former  case  corresponds  to  a trough  of  arbitreiry  cross 
section  and  has  already  been  considered.  The  latter  corresponds  to  a basin  with  rotational 
symmetry  where  Pq,  p„  are  now  the  polar  coordinates  of  source  and  receiver,  respectively, 
Tq,  measured  from  the  center  of  the  basin: 
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COS  e'n  sin  0;,  r„  = cos  Oq  sin  0o  '‘o-  (21) 

Thus  for  a circular  basin,  r cos  6 sin  0 is  invariant. 

For  a basin  of  any  other  shape,  there  are  two  alternatives:  use  Eq.  (19)  in  an  itera- 
tive fashion  since  the  integral  contains  dy  which  depends  on  ray  path,  or  divide  the  whole 
area  into  portions  of  circular  basins  or  strai^t  trouts,  each  with  its  own  separate, 
invariant,  matching  ray  angles  at  the  boundaries. 

Combining  Eq.  (12)  and  either  Eq.  (14)  or  (21),  we  can  find  a relation  for  azimuth 
angle  0„  with  a view  to  calculating  a continuous  curve  representing  the  horizontal  projec- 
tion of  the  ray  path.  After  many  surface  and  bottom  reflections,  iJ\e  difference  between 
®n»  ^'n  becomes  insignificant;  so,  dropping  the  subscripts  n and  primes,  we 

have  for  the  teough  case 

sin  0 = sin  0o  cos  Bq  HI(H^  - sin^  (22) 

The  distance  traveled  along  the  trou^  is  given  by 


X - j'dx  = J'dy  tan  0 = 


sin  00  cos  Bq  dy 


Hi 

(1  - sin^  00  cos^  Bq)  sin^  Bq 

i/2 


1/2 


(23) 


Given  the  cross  section  of  the  trough  H(y),  this  integral  can  be  performed,  and  the  ray 
path  y(ac)  is  then  known.  Some  examples  are  given  in  the  next  section.  Note  that  the 
trough  can  be  turned  into  a ridge  simply  by  making  i/(y)  an  increasing,  rather  than 
decreasing,  function  of  y. 

The  total  ray  length  is  given  by 


r.f.  ,rtan0dy  x 

s = I ds  = I dy  sec  0 sec  0 = I 2 — r^T"  = . 

J J J cos  Bq  sin  00  sin  0o  cos  Bq 

and  the  rather  surprising  result  is  that  the  ray  length  is  always  proportional  to  the 
component  of  distance  traveled  along  the  trough. 

For  the  circular  basin  (or  seamount),  Eq.  (22)  is  replaced  by 

sin  0 = sin  0o  cos  Bq  H y /(//2  - //§  sin2  0o)^^^. 


(24) 


(25) 
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and  the  ray  path  in  polar  coordinates  r,  4>  is  derived  from 


^ -j' tan  0 =j' 


sin  <t>Q  cos  Oq  dr 


' 2 f \ 

r — ( 1 ^ sin^  0Q ) - sin^  0o  cos^  Oq 

'•o  \ / 


The  total  ray  length  is  given  by 


s = J'ds  = J'dr  sec  0 sec  (p  = J 


tan  0 rdr J 

sin  (pQ  cos  0Q  Tq  Tq  sin  pQ  cos  0q 


2 (area  of  trajectory  swept  out) 

Tq  sin  pQ  cos  0Q 

The  ray  length  in  a basin  is  always  proportional  to  the  area  swept  out  by  the  radial 
coordinate  of  the  ray  as  it  moves  around  the  basin. 


II.  EXAMPLES  OF  HORIZONTAL  RAY  PATH  PROJECTIONS 
FOR  TROUGHS  AND  RIDGES 

A.  Trou^ 

Note  that  only  the  shape  of  the  tnugh  need  be  specified;  the  absolute  depth  is  not 
needed.  In  all  cases  the  characteristic  width  of  the  trough  is  R. 


0 fffyj  “ 11  +(y2/i?2)]-l/2 


''  (R^  + cot  0Q  sin  00 

(cot^  (?o  cos2  00  (R^  + y2)  + y2  _y2)l/2 


y = (cot^  0Q  cos^  00  {R^  + >„)  sin 


where  jc  is  a constant  given  by 


(R^  + cot  0Q  sin  0o 


[(R2  cot  00  sin  0oJ  [cot^  0o  cos^  pQ  {R^  + .Vo)  + ^0!^^^ 

So  the  rays  are  sine  waves  of  wavelength  2n{R^  + cot  0q  sin  Pq  and  amplitude 

[(R2  + y2)  cot2  00  COs2  pQ  + y§]l/2 
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Zero  depth 


Fig.  3 — Horizontal  projection  of  ray  paths  for  a 
trough  of  cosine  cross  section  H(y)  “ cos  (ylR), 
for  constant  tt>o,  and  9q  increasing  f Ti  zero  (indi- 
cated by  a broad  arrow) 


(ii)  H(y)  « COS  (y/R)  (cosine  trough) 


sin  (pQ  cot  Oq  cos  (y/R)  dy 


{[1  - sin^  00  cos^  Oq  '^0  ®o) 

(30) 


sin  (y/R)  = 1 


/ cos^  (yo//f)  sin^  00 

(x  x')(l  - sin^  00  cos^  0o) 

1 1 ~ 1 sin 

\ (1  - sin^  00  cos^  0o)/ 

R sin  00  cos  0o 

where  x'  is  given  by 

x'(l  - sin^  00  ^o) 


sin 


R sin  00  cos  00 


sin  (yo/R) 


cos^  (yo/R)  sin^  00 


(31) 


(1  - sin^  00  cos^  00 )/ 


When  the  first  term  on  the  right  of  Eq.  (30)  is  set  to  1 by,  for  instance,  making  yo  = 
nR/2,  which  corresponds  to  putting  the  source  at  the  edge  of  the  trough,  the  ray  path  is 
a periodic  zigzag  line  from  one  side  of  th.  trough  to  the  other  as  in  Fig.  3.  When  the 
same  term  is  made  small  compared  with  unity  by  making  yo/R  small  and  0o  or  <t>Q  large, 
the  ray  path  is  a sine  wave  of  wavelength 


27rft 


sin  00  cot  Oq 


(1  - sin^  00  cos^  Oq) 


and  amplitude 


/ cos^  (yo/^)  sin^  00 

R 1 , — 


(1  - sin'^  00  cos^^  Oq)/ 


L 
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(Hi)  H(y)  oc  /7  - (y^ /R^ )J  (elliptical  trough) 

(Note  that  the  small  bottom  slope  approximation  may  be  violated  near  the  edges.) 

r sin  00  cot  6^  (R^  - y2)l/2 

J (cot^  00  cos^  00  ■'■>'0 


X +x’  =Rsin  00  cot  Oq  E[k,  sin“^  (y/kR)] , (32) 

where 

x'  = R sin  00  cot  00  E[k,  sin”^  (y/feii)] , (33) 

and 

k = [cot^  00  cos^  00  + (y^//?2)]  1/2^  (34) 

where  E(k,  a)  is  an  elliptic  integral  of  the  second  kind 

E(k,a)=f  (1 -fe2  sin2  0)1/2 

•/o 

y is  still  always  periodic  in  x,  and  again  for  small  k;  that  is,  yo/R  small  and  0o  large  or  0o 
close  to  n/2,  E(k,  a)  = a,  and  the  ray  path  is  a sine  wave  of  wavelength  27rR  sin  0o  cot  0(; 
and  amplitude  Rk  = (R^  cot^  00  cos^  0o  + 


i 
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and,  for 


R>lyol> 


R cot  0Q  cos  00 
(1  + COt^  00  COS^  0o)^^^ 


y = (yo  - <^0*^  cos2  00  (R^  - y§))^/2  cosh  — — — (36) 

-yo)^  ^ cot  00  sin  00 

with  x'  given  by  putting  3C  = 0 with  y = yo-  These  curves  are  shown  in  Fig.  4.  The  cosh 
and  sinh  families  represent  rays  that  are  deflected  away  from  the  ridge  to  either  side. 

(ii)  H(y)  « cosh  (y/R) 

r sin  00  cot  00  cosh  (y/ii)  dy 


{[1  - sin^  00  cos^  00  - sin®  0o  cosh®  (yo/f^)]  + (1  - sin®  0o  cos®  0o)  sinh®  (y/R)}  t/® 

For  cosh  (Vq/R)  < (cosec®  0o  ~ sin®  0o  cot®  0o)^^^» 

/ cosh®  (yQ/R)  Y^®  (x  + x')(l  - sin®  0o  cos®  dn) 

sinh  (y/R)  = 1 - 1 sinh  : — 

V 0/1  • 0 • *0  r\  ’ ^ At. 

r 0Q  - sin*^  0Q  coi^  Oq 


yt/®  (x  + x')(l  - sin®  00  cos®  00 ) 

sinh  ;r— : — : ^ 

I R sm  00  cos  00 


and,  for  cosh  (yQ/R)  > (cosec®  0o  - sin®  0o  cot®  0o)^^^, 

/ cosh®(yo/f?)  \^^®  (jc  + jc')(l  - sin®  00  cos®  0o) 

cosh(y/fi)=  -1)  cosh  — — . 

y cosec®  6q  - sin®  0o  cot®  0q  j L ^ ^0  J 

(3 

These  curves  are  qualitatively  similar  to  those  of  the  last  example. 


Fig.  4 — Horizontal  projection  of  ray  paths  for  a 
ridge  of  cross  section  H(y)  ® (1  - 
for  constant  0q,  and  0g  increasing  from  zero  (in- 
dicated by  a broad  arrow) 
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Zero  depth 


Fig.  5 — Horizontal  projection  of 
ray  paths  for  tilted-plane  seabed 
ff(y)  “ y,  for  constant  (t>o,  and  Oq 
increasing  from  zero  (indicated  by  a 
broad  arrow) 


(Hi)  H(y)  oc  y (wedge) 


-/ 


sin  0Q  cot  Oq  y dy 


[(cosec-^  00  si"  ^o)  y -yol 


2 _«2il/2 


lx(\-  cos^  Oq  sin^  0o)  ^ 

y2  (1  - cos^  Oq  sin^  ^q)  = \ — 7 ^0  cos  Aq  cos  (pg 

\ sm  (pQ  cos  oq  ” " j 


+ yg  sin2  00 


The  ray  path  is  a hyperbola  as  shown  in  Fig.  5.  This  result  will  be  used  later  in  the 
intensity  estimation. 


III.  EXAMPLES  OF  HORIZONTAL  RAY  PATH 

PROJECTIONS  FOR  A BASIN  AND  A SEAMOUNT 

A.  Basin 

(i)  H(r)  oc  (1  + (r^/R^)}~^^^  (a  rotation  of  the  first  trough  cross  section) 

rgifi^  + cot  0Q  sin  (pg  dX 


*.r 

J 2X[-X^  + [(R2  +ri 


2X{-  X^  + [(R^  + r§)  cosec^  00  -R^]  Y - cot^  dg  sin^  0o 


where  X = and  fo  is  the  radial  coordinate  of  the  source. 
The  solution  of  this  integral  is 

A - B/r^ 


sin  (2(<h  + ‘h')l  = 


{A2  -2B)^'2 


I 

I 
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Fig.  6 — Horizontal  projection  of  ray  paths 
for  a circular  basin  of  cross  section  H(r  ) “ [ 1 + 
(r2/ii2) j-l/2^  forconstant0o. and  Oq  increas- 
ing from  zero  (indicated  by  a broad  arrow) 


where  4>'  is  a constant  given  by 


where 


sin  [2^'] 


A-B/rl 
{A^  - 


(41) 


A = -I-  Tq  ) cosec^  Oq  - 

B = 2 cot^  Oq  sin^  0o  ^q)- 

* 

The  ray  path  is  an  ellipse  centered  about  the  center  of  the  basin  (Fig.  6).  One  might 
have  expected  a rosette;  however,  for  this  particular  type  of  basin,  the  ray  covers  the 
same  path  again  after  one  revolution. 


An  ellipse  is  given  in  terms  of  its  minor  and  major  axes  a,  b,  respectively,  by 


m 
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The  ray  length  for  the  ellipse  (which  is  not  an  elliptic  integral  because  of  the  up-and- 
down  motion  of  the  ray)  is  given  in  terms  of  r by 


sin 


2rQ  sin  cos  Oq  (s  + s') 

(B/2)^/2 


2r2  ~A 

(i42-2B)^/2 


(42) 


or,  alternatively  in  terms  of  4>,  by 


tan 


Tq  sin  (pQ  cos  Oq  (s  + s") 

(B/2)1'2 


yj  tan  (4>  + <!»')  - (A^  - 2B)^^^ 


(43) 


where  s',  s"  are  given  by  substituting  s = 0 with  either  r = Tq  or  <I>  = 4>o.  For  a whole 
circuit  of  the  basin,  the  ray  length  is  given  from  Eq.  (42)  or  (43)  by 

Sc  = 27r(B^  + Tq)^/^  cosec  Bq.  (44) 

This  could  also  have  been  derived  by  substituting  wab  in  Eq.  (27)  for  the  area  of  ellipse 
swept  out. 


B.  Seamount 

(i)  H(r)  a r (conical  seamount:  apex  at  the  sea  surface) 


J r[r2  -r2(l 


r = rQ(l  - cos^  00  cos^  00 cosec 


(1  - cos^  00  cos^  00 


sin  00  cos  00 


(45) 


Here,  is  given  by  putting  4>  = 0 and  r = ro.  These  ray  paths  are  shown  in  Fig.  7. 


IV.  THE  EFFECT  OF  REFRACTION  ON  THE  RAY  PATHS 

Up  till  now,  we  have  assumed  the  sound  velocity  to  be  constant  so  that  surface-to- 
bottom  ray  segments  are  straight  lines.  This  is  probably  not  too  bad  an  approximation  in 
regions  where  there  is  surface  and  bottom  reflection  because  the  ray  angles  have  to  be 
comparatively  steep,  making  ray  curvature  in  the  vertical  plane  less  significant.  It  is 
possible  to  take  refraction  into  account  to  see  how  good  the  approximation  is. 


Equations  (1)  and  (4)  for  single  reflection  are  still  true  provided  the  angles  are 
measured  at  the  seabed.  Note  that  the  invariant 
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Sour 


Fig.  7 — Horizontal  projection  of  ray 
paths  for  a conical  seamount  H(r)  r, 
for  constant  ^o>  do  increasing  from 
zero  (indicated  by  a broad  arrow) 


derived  by  Weston  [1]  is  also  valid,  but  it  governs  the  integral  of  elevation  angle  over 
depth  as  opposed  to  the  elevation  angle  at  the  seabed.  The  first  alteration  comes  after 
Eq.  (9),  tvhere  the  spacing  of  adjacent  reflection  points  As,-  is  now  given  (with  z axis 
downwards)  by 


r^i  c(z)  cos  dj  dz 

As,- = 2 , (4 

“^0  (cf  - C^(z)  COS^ 

where  c,-  and  0,-  are  sound  velocity  and  incident  angle  of  the  ith  reflecting  point  on  the 
bottom. 


The  equivalent  of  Eq.  (6)  for  the  change  in  elevation  angle  at  the  /th  reflection 
point  is 


sin  0'i  = sin  0,-  / l - 


AHj  cot  0,- 
r c cos  0,  dz 

0 (cf  - cos^  0,)^ 


but  now  the  reflected  angle  at  the  (i  - l)th  point  is  related  to  the  incident  angle  at  the 
Ith  point  by  Snell’s  law: 


cos  0,-  = cos  0/_i  ( — ^ 
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The  total  change  in  angle  after  refraction  and  a single  reflection  is  then 

AHj  cot  dj 


sin2  0/  = 


f 


c cos  0,  dz 


(49) 


h (c?  - COS^ 


Aivjr  n rtilections,  we  have 


>2  " 


sin^  d'„  = sin^  0, 


1 - 


»(r)  n 

Vo/  /-I 

sCffMa 


AHi  cot  0i^ 


As, 


A«;COt0,Y 


As^. 


(50) 


Since  the  variations  in  sound  speed  c are  always  amall,  we  may  write  c = Cqq  + 5(z), 
and  Eq.  (46)  becomes 


. . r 

' C.  1 


H, 


As,  tan  0, 


(1  + 5/coo)dz 


' *^0  {1  - [(cqo  + 6)^  - c,?]  cot^  0,/c,?}  1/2 


(51) 


Now,  we  may  eliminate  0,-  from  the  correction  term  in  the  denominator  by  using  the  uni- 
form velocity  invariant  (Eq.  (12)): 


As,  tan  0,-  = 


2cr 


^00  p 


(1  + 6/coo)  dz 


'o  {1  - ((Cqo  + 5)^  - cf]  (cosec2  0(j  Hf/H^  - l)/c,?}  1/2 


(52) 


For  a given  velocity  profile,  As,  tan  0,-  is  simply  a function  of  Hq,  and  Oq,  and  so  we 
can  substitute  As,-  tan  0,-  into  Eq.  (50)  and  evaluate  the  products  as  we  did  in  the  earlier 
section  (assuming  that  2A//,-  cot  d,7As,-  is  always  small): 


sin2  0^  = sin2  0^  ( exp[-  2(F„  - Fq)] 

\%/ 


where 


t,  (\-^]exp[-2(F„-Fi)], 

i=i  \Ci  c._jy 


(53) 


H. 


Fn-Po 


/■ 


2dH 


As  tan  0 


1 


. I 


t i 
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The  summation  of  Eq.  (53)  may  be  evaluated  by  rearranging  its  terms  to  give 

Since  the  variations  in  c,-  are  always  small,  c,-  may  be  taken  out  of  this  summation  and 
Eq.  (53)  becomes 


sin2  0-  = gj„2  0^  ] e-2(^’n-^’o) 


J(l.a)-  P 

Vo 


where  a is  a number  very  close  to  unity  (of  order  c„/c,)  so  that  the  second  term  in  Eq. 
(54)  makes  only  a small  correction  and  the  uncertainty  in  a has  no  significant  effect. 
Equation  (54)  shows  that,  again,  the  elevation  angle  at  any  bottom  point  is  known  for  a 
given  initial  angle  and  water  depth.  Despite  the  more  complicated  nature  of  this  equation, 
contours  of  water  depth  on  a map  of  bottom  topography  are  still  contours  of  ray  eleva- 
tion angle. 

As  an  example,  a linear  velocity  profile  5 = [{z/Hq)  - 1]  Ac  gives  (assuming  ray 
angles  are  not  too  close  to  the  horizontal) 

r /i^.i  ie.iV 

A positive  velocity  gradient  (velocity  increasing  downwards)  results  in  a reduced 
spacing  of  reflection  points  for  a given  ray  elevation  angle  at  the  seabed,  so  a ray  travel- 
ing from  deep  to  shallow  water,  for  instance,  would  arrive  at  a given  elevation  angle 
sooner  (in  deeper  water)  than  it  would  if  it  had  traveled  through  water  of  uniform  velocity. 
As  expected,  the  refraction  effect  is  greatest  for  small  ray  elevation  angles. 

The  second  invariant  (Eq.  (14))  is  modified  only  slightly  because  each  reflection 
point  is  generally  at  a different  depth  so  that  the  elevation  angle  is  Eiltered  by  refraction. 
Incorporating  Snell’s  law  in  the  derivation  of  Eq.  (14),  we  have  for  a trough  of  arbitrary 
cross  section 

cos  e„  sin  <j>nlc„  = cos  Oq  sin  4>oIcq.  (56) 

This  equation  was  first  derived  by  Weston  [2] . Similarly,  for  a basin  with  rotational 
symmetry 


/I  • ^ 

— cos  sin  = — cos  Oq  sin  (pQ. 
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Clearly,  analytical  ray  paths  cannot  easily  be  obtained  using  Eqs.  (54)  and  (56)  or 
Eq.  (57),  but  numerical  solutions  are  possible. 


V.  SOUND  INTENSITY  IN  TROUGHS  AND  BASINS 

The  spreading  loss  at  a remote  receiver  may  be  calculated  numerically  by  finding 
the  final  spacing  in  three  dimensions  of  initially  adjacent  rays  as  given  in  sections  II  and 
III.  This  approach,  however,  does  not  throw  much  light  on  the  relative  importance  of 
loss  mechanisms  such  as  geometric  spreading  and  reflection  or  scattering  loss. 

A simpler  way  to  look  at  the  problem  of  propagation  loss  in  a basin  is  to  imagine 
(SOFAR)  rays  traveling  across  the  center  of  the  basin  at  constant  azimuth  because  of 
the  lack  of  bottom  interaction,  then  curving  round  in  a hyperbola  on  meeting  a locally 
sloping  plane  bottom  at  the  edge  of  the  basin  and  then  continuing  as  a SOFAR  ray  at  a 
new  heading.  If  we  take  the  wedge  model  of  Section  II  (ignoring  the  velocity  structure 
near  the  edge  of  the  basin),  we  see  that  the  asymptotes  of  the  hyperbola  intersect  at  the 
zero  depth  contour  (or  the  extrapolation  to  zero  depth  from  the  relevant  part  of  the  bot- 
tom if  the  bottom  slope  is  not  truly  constant),  so  the  horizontal  projection  consists,  to  a 
good  approximation,  of  straight  lines  reflecting  in  the  horizontal  plane  from  the  zero 
depth  contours  with  each  comer  rounded  off,  as  in  Fig.  8,  by  a hyperbola. 

Figure  9 shows  the  projection  of  two  slightly  diverging  sets  of  SOFAR  rays  arriving 
at  the  base  of  the  wedge  and  curving  round  in  almost  identical  hyperbolic  paths  so  that 
the  resulting  two  sets  of  SOFAR  rays  are  divergent  by  the  same  amount  as  before.  By 
the  invariants  derived  above  the  ray  elevation  angles  for  the  incoming  SOFAR  rays  on 
the  left  are  the  same  as  for  the  outgoing  rays  on  the  right.  Thus,  if  the  total  time  spent 
in  the  wedge  region  is  small  compared  with  that  spent  in  the  SOFAR  region,  the  part  of 
the  loss  associated  with  geometrical  spreading  is  the  same  as  if  the  SOFAR  ray  had 
carried  on  without  any  bottom  interaction  for  the  same  total  ray  length.  In  other  words, 
the  intensity  is  proportional  to  l/r  where  r is  the  total  horizontal  length  of  the  ray.  (It 
is  assumed  that  ranges  are  large  enough  for  the  effects  of  convergence  zones  to  have 
smoothed  out.) 


Fig.  8 — Method  of  constructing  the  hori- 
zontal projection  of  SOFAR  rays  after 
their  hitting  the  basin  edge 


V 
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In  practice,  of  course,  the  reflecting  surface  may  differ  from  a simple  tilted  plane, 
but  at  the  edge  of  a basin  the  greater  part  of  the  loss  would  probably  come  from  the 
individual  bottom  reflections  so  the  detailed  geometrical  effects  may  not  be  so  important. 


Zero  depth 


Fig.  9 — The  horizontal  divergence  d<l>  of  two 
sets  of  SOFAR  rays  is  practically  unchanged 
after  encountering  a straight  basin  edge 


The  number  of  bottom  reflections  encountered  near  the  edge  of  a trou^  is  given 
generally  by 


N 


=/- 


sin  9 ds  Ho  sin  00 


For  a wedge  of  slope  y. 


N 


_ _2  f 

“ >2  1 


2H  2 sin  <t)Q  cos  0o 


dy 


/tan  0 dy 
//2 


‘'yo  y [(cosec^  Bq  - cot^  0o  sin^  <PQ)y^  - y§]  ^2 
where  the  minimum  y coordinate  y„,  of  the  hyperbola  is 

= yo  sin  0o/(l  - cos^  6q  sin^  <t)Q)^'^. 

The  solution  of  the  integral  is 

sin  Oq 


^ ■ -1 
--sin 


(1  -cos^  00  sin^ 


(58) 


(59) 


(60) 


where  the  fact  that  y = H^lyQ  has  been  used.  Provided  that  6q  is  small  and  <1>q  is  not  too 
close  to  n/2,  this  may  be  expanded  as 


N=  i -00  sec0o)-  (61) 

Since  most  SOFAR  rays  are  likely  to  hit  the  sloping  bottom  when  they  are  at  the 
lowest  point  of  their  trajectory,  their  initial  ray  angles  for  the  above  calculation  are 
likely  to  be  close  to  zero  despite  the  fact  that  the  initial  angles  for  the  SOFAR  rays  were 
arbitrary.  The  order  of  magnitude  of  loss  L after  an  encounter  with  the  edge  of  the  basin 
is  given  in  terms  of  R,  the  loss  per  bottom  reflection  (a  logarithmic  quantity)  by 
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L = NR=  ^ R. 

27 

R is  in  reality  different  for  each  reflection  since  it  is  a function  of  elevation  angle;  how- 
ever, let  us  assume  the  worst  case  (the  steepest  angle)  where,  since  at  this  point  (the 
closest  approach  to  shallow  water)  sin  ^>  = 1,  cos  Qm  ~ cos  9q  sin  (pQ.  As  we  have  already 
assumed  6q  to  be  small,  we  have  9^  ^ 7r/2  - <I>q. 

Thus,  the  appropriate  grazing  angle  in  the  vertical  plane  for  evaluation  of  the  reflec- 
tion coefficient  is  simply  the  horizontal  grazing  angle  to  the  edge  of  the  basin  (Fig.  8) 

R=R(9)==^R(v/2-(t>o).  (63) 

The  total  loss  on  encountering  the  basin  edge  may  be  thought  of  as  the  loss  given  by 
a single  specular  reflection  at  a vertical  surface  along  the  zero  depth  contour*  made  of  the 
same  material  as  the  bottom  multiplied  by  the  real  number  of  reflections  (which  is  roughly 
n/2  divided  by  the  bottom  slope).  This  is,  of  course,  equivalent  to  the  empirical  scattering 
cross  section  already  mentioned.  The  same  interpretation  serves  as  a quick  method  of 
constructing  horizontal  ray  paths,  and  it  may  be  shown  by  use  of  Eq.  (24)  that  ray  lengths 
are  correct  as  well. 

As  an  example,  the  reflection  loss  per  bounce  near  the  edge  of  the  basin  with  = 
70°  as  in  Fig.  8 (and  consequently  9^  = 20°)  could  be  1 dB  with  the  20°  grazing  angle 
and,  say,  a source  frequency  of  10  Hz.  A bottom  slope  of  2°  would  give  45  bottom 
reflections  and  an  overall  reflection  loss  of  45  dB.  The  only  other  loss  is  that  due  to  the 
usual  spreading  law  for  a source  and  receiver  spaced  apart  by  the  total  real  ray  length. 

If,  instead  of  the  tilted-plane  reflecting  surface,  a surface  with  gently  curving  depth 
contours  such  as  a seamount  had  been  chosen,  the  main  difference  in  the  above  calcula- 
tions would  have  been  in  the  horizontal  divergence  of  outgoing  rays.  From  simple  two- 
dimensional  geometry  it  can  be  shown  that  the  intensity  vs  range  relation  is  now 


ICC 


u 

ri(r2  + u) 


where 


r j p cos 

p cos  - 2rj  ’ 


(64) 


(65) 


where  rj,  r2,  and  p are  source-seamount  and  seamount-receiver  ranges  and  seamount 
radius  of  curvature,  respectively,  and  2<P  is  the  overall  change  in  bearing  of  the  ray.  As 
expected,  this  makes  the  intensity  inversely  proportional  to  the  total  range  (rj  + r2)  as 
the  seamount  radius  of  curvature  tends  to  infinity,  and,  as  the  curvature  becomes  small 
compared  with  rj  or  r2,  the  intensity  is  proportional  to  (p/2)  cos  4>/(r^r2)  which  is 
equivalent  to  Weston’s  target  strength  formula. 


♦In  the  more  general  case  of  a bottom  with  unequally  spaced  contours,  this  would  have  to  be  the  zero 
depth  contour  at  the  tangent  plane  to  the  bottom  at,  say,  the  closest  approach  to  shallow  water. 
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VI.  CONCLUSION 

Ray  invariants  were  derived  and  used  to  calculate  analytically  the  horizontal  projec-  j 

tion  of  surface-  and  bottom-reflected  ray  paths  in  the  vicinity  of  submarine  troughs,  I 

ridges,  basins,  and  seamounts  of  various  cross  sections.  The  types  of  trajectories  derived 

include  a sine  jwave  in  a trough,  a hyperbolic  sine  and  hyperbolic  cosine  at  a ridge,  a ; 

hyperbola  at  a wedge,  an  ellipse  in  a circular  basin,  and  a hyperbola-like  curve  at  a conical 
t seamount. 

The  ray  length  in  a trough  of  constant  cross  section  was  found  always  to  be  propor- 
tional to  the  component  of  distance  traveled  along  the  trough,  and,  in  a basin  with  rota- 
tional symmetry,  the  ray  length  was  proportional  to  the  area  swept  out  by  the  radial 

coordinate  of  the  ray.  ; 

Refraction  was  found  not  to  affect  greatly  the  horizontal  projection  of  the  ray  path, 
although  it  makes  solution  in  terms  of  simple  functions  impossible. 

The  reflection  loss  and  geometrical  spreading  at  the  edge  of  a basin  behave  as  if  | 

I there  had  been  a vertical  reflecting  surface  at  the  zero-depth  contour  (see  footnote,  p.  19),  ' 

with  reflection  coefficient  identical  to  that  of  the  real  bottom  except  for  being  raised  to  I 

the  power  7r/(2  bottom  slope).  j 

Many  authors  have  adequately  modeled  echoes  from  seamounts  and  basin  edges  with  ! 

an  empirical  scattering  cross  section  or  target  strength  for  each  feature.  Some  authors  j 

j include  a detailed  single-scattering  calculation  [6] ; others  make  no  attempt  to  explain  the  ! 

target  strength  in  detail  and  only  imply  a scattering  (as  opposed  to  reflection)  model  ^ 

[7-9].  Only  very  rarely  are  the  data  entirely  incompatible  with  the  type  of  multiply- 
reflected  ray  paths  investigated  in  this  report,  and  in  most  cases  these  ray  paths  would 

probably  fit  the  data  more  closely.  In  practice,  determination  of  the  correct  mechanism  i 

by  measuring  the  propagation  loss  is  likely  to  be  hampered  by  poor  knowledge  of  i 

environmental  data,  such  as  bottom  roughness  and  reflection  coefficient.  Despite  this, 
an  experiment  could  be  designed  specifically  to  distinguish  between  the  two  mechanisms 
by  measuring  different  variables  such  as  transmission  angle  and  arrival  angle  in  the  hori-  i 

zontal  plane.  j 
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